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An  important  task  of  the  Optimum  Pumping  Technology  program  is  the 
investigation  of  the  feasibility  of  discharge  pumping  electronic 
transition  lasers.  Such  an  investigation  requires  a detailed 
knowledge  of  the  secondary  electron  transport  properties  that  may 
be  obtained  from  a numerical  solution  of  the  Boltzmann  code.  An 
example  of  an  electronic  transition  discharge  pumped  laser  is  HgCl* 
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20.  Abstract  (continued) 

A theoretical  investigation  of  the  e-beam  controled  discharge  pump- 
ing of  the  mercury  halide  lasers  requires  a coupled  solution  of  th  5 
electron  transport  Boltzmann  Equation  and  dominant  kinetics  rate 
equations.  The  general  only  the  lowest  order  diffusion  approxima- 
tion of  the  Boltzmann  Equation  is  used  for  this  purpose.  This  is 
because  the  more  general  Boltzmann  Equation  is  complex  in  nature. 

To  insure  the  validity  of  the  approximation  an  analysis  to 
determine  the  error  of  the  lowest  order  diffusion  equation  is  made 
The  analysis  shows  that  if  the  spatial  derivative  of  the  distribu- 
tion function  can  be  neglected,  i.e.,  for  a distributed  source  of 
electrons,  the  Pi  approximation  is  valid  for  all  the  cases  con- 
sidered. When  the  spatial  derivative  cannot  be  neglected,  which 
is  the  case  for  a spatially  localized  source  of  electrons,  then  the 
region  of  validity  of  the  Pi  approximation  is  much  more  restric- 
tive. High  pressure  discharges  are  examples  of  a distributed 
source.  Electron  beams  generated  for  swarm  experiments  are  ex- 
amples of  a spatially  localized  source. 

In  order  to  address  the  key  technical  issue  of  discharge  stability 
in  mercury  halide  lasers;  a theoretical  study  is  made  to  predict 
the  stability  criteria  in  recombination  and  attachment-dominated 
metal  vapor  discharges.  In  mercury  halide  laser  discharges,  a 
Relatively  large  excited  state  population  exists  which  dominates 
the  stability  of  these  discharges.  This  occurs  because  the  ex- 
cited states  have  a smaller  ionization  potential  and  larger  ioniza- 
tion cross-section  than  the  ground  state.  As  a result,  multi-step 
ionization  is  the  dominant  ionization.  When  three  step-ionization 
is  dominant  and  electron  mixing  of  the  excited  states  is  unimpor- 
tant, attachment-dominated  discharges  will  be  stable  if  the  attacl 
ment  rate  8>m  vig  where  2<m<3  and  v^g  is  the  equilibrium 
ionization  rate.  For  a recombination  dominated  discharge  the  cri- 
terion in  this  limit  is  a n6g  > m vig  where  1.0<m<1.5  and  a is 
the  recombination  coefficient  and  n^n  is  the  equilibrium  electron 
density.  0 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfNTlcn  Dtlm  Enfrtd) 


TABLE  OF  CONTENTS 


Section  Page 

I . INTRODUCTION  3 

II.  ON  THE  DIFFUSION  APPROXIMATION  OF  THE  BOLTZMANN  7 

EQUATION  FOR  SLOW  ELECTRONS 

III.  STABILITY  CRITERIA  IN  METAL  VAPOR  DISCHARGES  19 

REFERENCES  31 


1 


sdOKJ CO  EVERETT 

\J  '-J 


In  the  last  three  years  ARPA/ONR  sponsored  Visible/UV  laser 
programs  of  AERL  have  led  to  the  discovery  and  successful  devlop- 
ment  of  important  rare  gas  halide  lasers.  Specifically,  we  have 
obtained  12  j/liter  with  a 10%  intrinsic  efficiency  from  KrF*  by 
pure  e-beam  pumping  and  10  J/liter  with  a 9.5%  intrinsic  efficien- 
cy from  e-beam  controlled  discharge  pumping.  In  the  last  year 

another  class  of  visible  lasers,  the  mercury  halides,  have  been 
(1  2 ) 

discovered.  ' These  molecules,  like  the  rare  gas  halides,  have 
ionic  upper  levels.  Hence,  the  formation  kinetics  of  the  upper 
laser  level  should  be  rapid  and  efficient  due  to  the  long  range 
attractive  forces  involved  in  ion-ion  interactions.  Since  the 
ionization  energy  of  mercury  is  lower  than  that  of  the  rare  gases, 
the  mercury  halides  radiate  at  a longer  wavelength  than  the  rare 
gas  halides.  To  date,  lasing  action  has  been  obtained  in  HgCl 
at  5576  8 and  HgBr  at  5018  8. 

E-beam  controlled  discharge  pumping  of  the  mercury  halide 
lasers  could  lead  to  very  high  overall  laser  system  ef f iciencies- 
efficiencies  as  high  as  20%.  In  order  to  achieve  these  high 
efficiencies,  efficient  production  of  mercury  3p  manifold, 

Hg*  ("*P)  , in  a stable  discharge  run  at  high  enhancement  ratios 
(discharge  power/e-beam  power)  must  be  demonstrated.  The  key 
technical  issues  are: 
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1)  Hg  (^P)  production  efficiency 

2)  Discharge  stability 

3)  Discharge  enhancement  ratio 

In  order  to  theoretically  predict  the  characteristics  of 
the  mercury  halide  lasers  a coupled  solution  of  the  electron 
transport  equation  and  kinetic  rate  equations  that  model  the  key 
formation  and  quenching  kinetics  is  required.  Such  a solution 
is  necessary  because  the  electron  impact  excitation  and  ioni- 
zatoin  rates  are  dependent  on  the  discharge  power  deposited  in 
the  gas.  The  solution  of  the  Boltzmann  Equation  for  this  purpose 
is  generally  obtained  by  considering  only  the  lowest  order  dif- 
fusion or  approximation.  This  is  because  the  more  general 
form  of  the  Boltzmann  Equation  is  not  amenable  to  analytic  solu- 
tions. The  P^  approximation  of  the  Boltzmann  Equation  has  proved 
most  useful  in  the  understanding  of  high  pressure  discharges. 
However  such  an  approximation  is  not  necessarily  valid  for  all 
cases  under  consideration.  Hence,  the  first  part  of  this  report 
investigates  the  regions  of  validity  of  the  P^  approximation. 

One  of  the  major  technical  issues  of  the  e-beam  controlled 
discharge  pumping  of  the  mercury  halide  lasers  is  the  discharge 
stability.  The  mercury  halide  laser  discharges  which  are  typi- 
cal of  attachment-dominated  discharges  are  characterized  by  a 
relatively  large  excited  state  population.  These  excited  states 
dominate  the  stability  of  these  discharges.  Typically  they  have 
a smaller  ionization  potential  and  a larger  ionization  cross 
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section  than  the  ground  state.  As  a consequence  these  effects 
lead  to  a condition  where  multi-step  ionization  dominates.  A 


theoretical  study  to  establish  the  stability  criteria 
attachment  and  recombination  dominated  discharges)  is 
in  the  second  part  of  this  report. 
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II.  ON  THE  DIFFUSION  APPROXIMATION  OF  THE 
BOLTZMANN  EQUATION  FOR  SLOW  ELECTRONS 

The  lowest  order  diffusion  or  P^  approximation  of  the 

Boltzmann  equation  has  proved  most  useful  in  the  understanding 

of  high  pressure  discharges.  There  have  been  several  efforts 

addressing  the  domain  of  applicability  of  this  approximation.^ 

( 2 ) 

In  one  of  the  earliest  discussion  Holstein  has  shown  that  the 
diffusion  approximation  is  valid  when  the  elastic  cross  section 
for  electron  heavy  particle  collisions  is  large  compared  to  the 
inelastic  cross  section.  This  criterion  has  also  been  used  in 
the  Bethe  Age  theory  for  electron  and  neutron  diffusion. ^ ^ While, 
for  most  gases,  this  criterion  is  well  satisfied,  there  are  in- 
stances when  the  magnitude  of  inelastic  and  elastic  cross  sections 

( 5 ) 

are  comparable.  For  example,  recent  measurements  of  the  vibra- 
tional cross  section  of  N^  evaluates  the  peak  cross  section  in  the 
energy  range  of  2-3  eV,  to  about  12  R,  which  is  comparable  to  the 
peak  elastic  cross  section  of  30  i?2  for  N2  evaluated  by  Englehardt, 
Phelps  and  Risk.  ^ 

In  this  article  we  investigate  the  applicability  of  the  P^ 
approximation  by  considering  a narrow  band  of  electrons  in  energy 
space.  It  is  shown  that  the  angular  distribution  function  depends 
on  two  dimensionless  parameters,  y and  y is  the  ratio  of  the 

energy  gained  by  an  electron  in  the  electric  field  per  transport 
mean  free  path  to  the  electron  energy  and  £ is  the  ratio  of  the 
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energy  lost  by  inelastic  collisions  per  transport  mean  free  path 
along  its  trajectory  to  the  electron  energy.  The  former  is  a 
familiar  parameter  in  plasma  physics,  while  the  later  is  propor- 
tional to  the  ratio  of  the  inelastic  to  the  elastic  cross  sections. 
The  analysis  presented  suggests  that  when  the  spatial  derivative 
of  the  distribtuion  function  cannot  be  neglected,  the  approxi- 
mation is  not  valid  for  values  of  y and  E,  that  approach  unity. 
However,  when  the  spatial  derivative  can  be  neglected,  the  error 
in  Pj  approximation  is  limited  to  5%  for  all  electron  energies 
and  for  all  values  of  y and  E,  considered.  These  results  imply 
that  the  applicability  of  the  P^  approximation  depends  strongly 
on  whether  the  sources  of  electrons  are  spatially  localized  or 
uniformly  distributed  throughout  the  volume  under  consideration. 
This  is  because  for  the  localized  source,  the  spatial  variation 
of  the  electron  density  is  important,  while  for  a uniformly  dis- 
tributed source  the  spatial  variation  may  be  neglected.  Spatially 
localized  sources  of  electrons  that  diffuse  between  two  electrodes 
are  used  in  swarm  experiments.^  Hence,  for  gases  like  for 

which  values  of  y and  E,  are  large,  the  evaluation  of  impact  param- 
eters from  swarm  data  using  the  P^  appoximation  will  be  suspect. 

In  a high  pressure  discharge,  however,  secondary  electrons  are 
created  uniformly  (at  least  on  the  scale  length  of  the  transport 
mean  free  path)  through  the  discharge  volume.  For  such  cases, 
therefore,  the  spatial  derivative  of  the  distribution  function 
may  be  neglected. 
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The  transport  of  a broad  area  electron  source  in  the  presence 


of  collisions  can  be  described  by  the  one  dimensional  Boltzmann 

(8) 


equation , 


0 , 3f  0 „ , 3f  2.  3f 

2y  A ^ + 2 eE  A U ^ + Y (1  — M ) ^ = 2A 


3f 


3s 


(1) 


where  the  distribution  function  f is  a function  of  the  kinetic 
energy  w,  angular  coordinate  y and  spatial  coordinate  x.  p is 
the  cosine  of  the  angle  between  electron  trajectory  and  E,  which 
is  assumed  to  be  along  the  x direction.  The  transport  mean  free 
path  A is  defined  in  the  usual  manner.  The  parameter  y is  the 
ratio  of  the  energy  gained  per  transport  mean  free  path  to  the 
electron  energy,  i.e., 

eE  A 


Y = 


w 


(2) 


3 f 


For  isotropic  elastic  scattering,  -r— 

o S 


is  given  by 


(8) 


+ i 


= -2f  + f f ( y 1 ) dy 
e -1 


(3) 


We  will  assume  that  an  inelastic  scattering  process  changes  only 


th 


the  energy  of  the  electron.  Then  for  the  J excitation  process 


3f 
3 s 


Pj  (w  + Ej)  f (w  + Ej) 


Pj  (w)  f (w)  | 


(4) 


where  Pj (w)  is  the  average  number  of  collisions  per  unit  distance 


along  the  electron  trajectory.  E T is  the  energy  lost  after  the 

J 


collision.  PT(w)  is  zero  for  w < E . 

J J 
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The  usual  method  of  solving  the  Boltzmann  equation  for 


high  pressure  discharges  is  to  neglect  the  variations  of  f with 
x and  then  expand  f as  a sum  of  Legendre  polynomials.  Because 
of  the  complexity  only  the  first  two  polynomials  are  generally 
retained.  This  procedure  is  valid  provided,  of  course,  that  f 
is  nearly  isotropic.  A method  for  checking  the  validity  of  such 
an  approach  is  to  simplify  Eq.  (1)  by  considering  a narrow  beam 
of  electrons  (of  width  Aw)  in  energy  space  having  a mean  energy 
Wg.  Integrating  Eq.  (1)  over  all  energy,  we  obtain 

+1 

2uAo  H + 2v  - ^ovF  + v1  - »2>  If  * - 2F  + f f«>*»  a“ 


where 


oo 

-f 


dw  w f (w) 


The  dimensionless  parameter  £ is  the  ratio  of  energy  lost  by  the 

electrons  per  transport  mean  free  path  to  the  energy  of  the  elec- 
(8) 


tron 


= when  EJ 


Aw  <<  w. 


^0^0 

C0  = -E-  when  Ej  ^ w0  (8) 

The  two  definitions  of  are  obtained  from  Eq.  (6)  for  the  two 
limiting  conditions.  Equation  (7)  is  obtained  by  a Taylor  expan- 
sion of  the  first  term  on  RHS  of  Eq.  (4) , while  for  a narrow  beam 
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■ 


this  term  will  have  a negligible  contribution  when  « w.  Xg 
is  the  energy  lost  by  the  electron  per  unit  path  length  along  its 
trajectory  and  is  given  as 


‘o'E 


EJPJ  = N 


QJEJ 


(9) 


.th 


where  Q is  the  inelastic  cross  section  for  the  J process  and 

vJ 

N is  the  number  density  of  scatterers.  From  Eqs.  (7) , (8)  and 

(9) , it  is  clear  that  £q  is  proportional  to  the  ratio  of  the 
inelastic  cross  section  to  the  elastic  cross  section.  In  fact, 
when  Ej  ~ Wq , which  is  near  or  just  above  the  threshold  of  the 
excitation  process,  is  precisely  the  ratio  of  the  inelastic 
cross  section  to  the  elastic  cross  section.  Note  that  eval- 
uated by  Eq.  (7)  is  always  less  than  when  evaluated  by  Eq.  (8) , 
except  when  wQ  = Ej.  As  the  two  values  are  similar  and  we  are 
only  making  estimates  of  the  error,  we  will  use  Eq.  (7)  in  the 
subsequent  discussion.  y^ , and  Xg  are  all  evaluated  at  the 

particular  energy  w^. 

Dropping  t’..e  subscript  0 and  integrating  Eq.  (5)  over  all 
y we  get 


, dJ  FJ_ 

A dx  + <y>  “ u 


(10) 


where  J is  the  electron  current  density  and 


+ 1 

r ) 

<y>  = 

/ Fydy  J 

/ Fdy  J 

J-1  l -1  ' 
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— 1 


J = JQ  exp  - (/  dx’  ^ A, 


Assuming 


F ~ g (y ) exp  (-  / dx ' A) 


Eq.  (5)  reduces  to 


+ 1 

+ 2C  g - 2y  yg  + y (1  - y2)  = _ 2g  + J g(y')  dy ' 


Making  a Legendre  polynomial  expansion  of  g(y) 


g (y ) = Z A£P£ (y)  (14) 

Equation  (13)  can  be  represented  by ’ an  infinite  set  of  simultaneous 
equations 

Vl  + \ + 2«1 

- An-1  <2^'  + <n+1>  If]*  0 (15) 

Using  Eq.  (15)  we  may  solve  for  the  coefficients  A£  by 
truncating  the  series  at  some  arbitarary  SL  = L.  An  iterative 
procedure  is  necessary  as  <y>  is  not  known  a priori.  In  this 
manner  one  may  obtain  a solution  of  g(y)  to  any  desired  accuracy. 
Note  that  when  the  variation  of  f with  x is  neglected  then  this 
equation  takes  a simpler  form  (n  > 0) 
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for  which  an  iterative  procedure  is  not  necessary.  The  error  in 
the  approximation  may  be  evaluated  by  the  following  ratio 

An2/  ^ An2  (17) 

n = 2 n = 0 

It  is  interesting  to  note  that  the  series  represented  by  Eq.  (16) 
converges  much  more  rapidly  than  Eq.  (15)  for  a given  y and  £, 
signifying  that  when  the  spatial  derivative  is  neglected  the  dis- 
tribution function  is  much  more  isotropic.  The  error  for  the 
case  when  electron  distribution  is  not  a function  of  x,  is  < 5% 
for  all  y and  C values  considered  in  Figure  1.  However,  the  same 
argument  does  not  apply  to  the  situation  when  the  variation  of  f 
with  x cannot  be  neglected,  as  in  the  case  of  electron  swarm  ex- 
periments. This  result  is  supported  through  Figure  1 that  shows 
the  contours  of  the  error  as  a function  of  y and  £.  As  one  ex- 
pects, the  error  increases  as  both  y and  £ increase.  Figure  2 
shows  the  distribution  function  obtained  from  our  Boltzmann  code 
along  with  our  estimate  of  error  and  the  parameters  y and  £ for 
an  electric  field  of  1 kV/cm  amagat  in  nitrogen.  The  error  < 1% 
for  most  of  the  electrons  having  an  energy  > 0.2  eV.  The  error 
is  large  for  a few  electrons  near  zero  energy  as  Holstein  has 
pointed  out.  This  is  because  y is  large  in  this  region.  Since 
relatively  large  numbers  of  electrons  have  energies  beyond  0.2  eV 
for  this  electric  field  (as  is  observed  from  the  distribution 
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Figure  1 % Error  in  Using  Approximation  for  Various  Values 

of  Parameters  y and  £ 
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Figure  2 % Error  in  Using  Approximation  for  Nitrogen  at 

1 kV/cm  Amagat  at  Various  Electron  Energies 


15 


^27AVCO  EVERETT 


function) , the  approximation  would  be  a valid  approximation 
for  this  case.  However,  when  the  electric  field  is  increased  to 
20  kV/cm  amagat,  the  value  of  y increases  by  a factor  of  20  and 
the  corresponding  error  and  calculated  value  of  E,  is  shown  in 
Figure  3.  It  is  observed  that  the  error  in  using  approxima- 
tions for  this  case  is  quite  high. 

The  results  in  this  article  indicate  that  the  evaluation 
of  electron  impact  parameters  from  swarm  data  using  the  P^  ap- 
proximation will  be  suspect  in  gases  like  N2*  However,  if  the 
detailed  cross  sections  are  obtained  from  beam  experiments,  these 
include  the  momentum  transfer  cross  section,  the  P^  approximation 
is  valid  for  high  pressure  discharges. 
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Figure  3 


% Error  in  Using  Approximation  for  nitrogen  at 
20  kV/cm  Amagat  at  Various  Electron  Energies 
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III.  STABILITY  CRITERIA  IN  METAL  VAPOR  DISCHARGES 

The  stability  criterion  for  discharges  where  excited  state 
(two-step)  ionization  plays  rn  important  role  has  been  previously 
investigated.^  Such  a stability  criterion  has  been  shown^*^ 
to  describe  high  power  rare  gas  halide  laser  discharges  when  the 
discharge  electric  field  is  held  constant  during  the  discharge 
pulse.  In  this  letter,  the  stability  of  metal  vapor  discharges 
is  investigated.  Interest  in  metal  vapor  lasers  has  been  height- 
ened since  the  observation  of  lasing  action  in  HgCl*  at  5576  £ 

by  both  pure  e-beam  pumping ^ ^ ^ and  e-beam  controlled  discharge 
(12) 

pumping.  In  the  high  power  metal  vapor  discharges,  there 

• ( 1 2 ) 

exists  a relatively  large  excited  state  population  which 

dominates  the  stability  of  these  discharges.  This  occurs  becau  e 

typically  the  excited  states  have  a smaller  ionization  potential 

and  a larger  ionization  cross  section  than  the  ground  state. 

For  high  power  discharges  such  as  those  used  in  pumping  electronic 

transition  lasers,  these  effects  lead  to  a condition  where  multi- 

step  ionization  dominates.  Both  attachment  and  recombination 

metal  vapor  discharges  have  been  considered  in  this  paper.  The 

mercury  halide  lasers  are  examples  of  attachment  dominated  metal 

(14) 

vapor  discharges.  The  CdHg*  laser  is  typical  of  recombination 

dominated  metal  vapor  discharges.  In  the  subsequent  discussion 
on  discharge  stability,  conditions  typical  of  those  found  in  these 
lasers  will  be  considered  as  examples  of  the  stability  criteria 
derived. 
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Figure  4 shows  the  energy  levels  of  cadmium,  mercury  and 

argon  (used  as  a buffer  gas).  In  cadmium  for  example,  it  requires 

3.7  eV  to  create  Cd*  and  another  5.3  eV  to  ionize  Cd*  whereas  it 

requires  only  2.6  eV  to  produce  Cd**  from  Cd*  or  to  ionize  Cd**. 

As  indicated  in  Figure  4,  Cd*  and  Cd**  are  a number  of  relatively 

closely  spaced  levels.  As  the  electron  temperature  will  typically 

be  larger  than  this  spacing,  we  assume  that  the  manifold  of  states 

represented  by  Cd*  and  Cd**  will  be  completely  mixed  and  hence 

can  be  represented  by  a single  level.  Since  the  average  electron 

energy  in  these  discharges  is  1-1.5  eV  three-step  ionization 

tends  to  dominate  at  the  electron  densities  of  interest.  However, 

if  the  electron  density  is  large  such  that  the  mixing  of  the  two 

excited  states  (Cd*  and  Cd**)  by  collisions  of  the  first  and 

second  kind  is  far  more  rapid  than  excited  state  reactions  with 

other  heavy  particles  or  ionization  then  the  Cd*  and  Cd**  may  be 

(9) 

treated  as  a single  level  and  the  two-step  ionization  criterion 
derived  previously  is  valid.  In  the  cadmium  mercury  system, 
stabilization  of  the  discharge  current  will  be  achieved  through 
loss  of  secondary  electrons  by  recombination.  In  the  case  of  the 
HgCl*  laser,  the  excited  states  Hg*  and  Hg**  will  play  a dominant 
role  in  discharge  stability  and  discharge  stabilization  will  be 
achieved  through  loss  of  electrons  by  attachment. 

Three  equations  are  important  in  determining  the  stability 
of  the  discharge  in  the  presence  of  rapid  excited  state  ionization. 
The  first  describes  the  production  and  loss  of  discharge  electrons. 
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Figure 


CADMIUM  MERCURY  ARGON 

4 Energy  Levels  for  Cadmium,  Mercury  and  Argon 
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dF  = seb  + ,vi  - S)ne  «18> 

while  the  second  describes  the  production  and  loss  of  HgMCd*), 
dn.. 

dF  = K*  nen  ' V*  nenl  + V*  nen2  " Ranl  " V nlne  <19) 

and  the  third  describes  the  production  and  loss  of  Hg**(Cd**)f 
dn 

dF  " V*  nenl  - V*  nen2  * Rbn2  " vine  <20> 

Assuming  that  the  dominant  ionization  is  the  ionization  of 

v.  is  defined  as 
1 

v±  = K+  n2  (21) 

where  K+  is  the  ionization  rate  constant  of  n2  by  electron  impact. 
Here  ng  is  the  electron  density;  Seb  is  the  electron  production 
rate  resulting  from  the  high  energy  e-beam;  3 is  the  attachment 
rate.  Note  that  3 = ang  for  recombination  dominated  laser  (CdHg*); 
a is  the  recombination  rate;  K*  is  the  production  rate  constant  of 
n^^  by  electron  impact;  K^**  is  the  production  rate  constant  of  n2 
by  electron  impact;  K^*  and  K^**  are  the  superelastics  rate  con- 
stant of  n^  and  n2  by  electron  impact  respectively;  Rq  is  the 
inverse  lifetime  of  n^;  is  the  inverse  lifetime  of  n2;  n,  n.^ 
and  n2  represent  the  number  densities  Hg(Cd),  Hg*(Cd*)  and 
Hg**(Cd**)  respectively.  Figure  4 depicts  schematically  the 
terminology  of  the  rates  described  above.  In  the  present  ana- 
lysis we  will  assume  that  the  ground  state  density  n is  constant. 
Further,  only  transitions  between  adjacent  levels  have  been  con- 
sidered as  these  rate  constants  will  be  by  far  the  most  rapid. 
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Equations  (18)  to  (20)  are  three  nonlinear  simultaneous 
differential  equations.  If  the  discharge  electric  field  is  a 


constant  in  time,  we  can  linearize  by  letting, 
ne  = nQ  + ne(t),  n1  = n1  + n1(t),  n2  = 


n = n + n (t) , n,  = n + n,  (t)  , n9  = n_  + n9(t)  (22) 

e eQ  e x iQ  x z zQ  z 


The  equilibrium  conditions,  n , n1  and  n„  is  readily  obtained 

e0  x0  0 

from  Eqs.  (18)  to  (20).  The  normalized  perturbed  Eqs.  (18)  to 
(20)  are  given  by, 

dn  _ _ 

~dt"  = -Clne  + C2n2  (23‘ 


dn^ 

dt 

= R n + 

C,n, 

a e 

3 . 

dn2 

ll 

3 

+ 

Ccn 

dt 

b e 

5 

= n-^/nj^  , 



ni 

n2 

an  + S , 
e eb 

C,  = Recombination  dominated 

1 n 

e0 


Attachment  dominated 


(26a) 


(26b) 


Kd**  ne  n2 

C2  + K+  nv  C3  = —2 0,  C4  = Ra  + Ku*‘  n^,  Cj 

' Kd**  neQ  + Rb 


(27a-d) 


Equations  (23),  (24)  and  (25)  can  be  combined  by  elimination  to 

yield 
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d3  n 


d2  n 


— + B -t—  + Cn  =0 
dt  e 


where  A,  B and  C are  functions  of  C^,  C 3,  and  C For 

Eq.  (28)  to  be  of  stable  type  (i.e.,  the  characteristic  equation 
has  real  negative  roots  or  the  real  part  of  the  complex  roots  be 


negative)  the  following  conditions'  must  be  satisfied, 

A>0,  B>0,C>0,  AB-C  > 0 


(29a-d) 


We  have  analyzed  Eq.  (29a-d)  in  order  to  establish  stability 
criteria  for  recombination  and  attachment  dominated  discharges. 

The  most  restrictive  condition  is  obtained  from  Eq.  (29c)  (i.e.. 


C > 0)  and  is  given  as. 


an  > (1  + e/2)  v.  Recombination  dominated 
e0  10 


(30a) 


B > (2  + e)  v. 

1 / 


Attachment  dominated 


(30b) 


+ A6162_1 
where  v.  = K n , e = - -t-  . , 

iQ  2q  l+A(61+62)+A6162 


61  (K  **  + K *)n 
u d e. 


r _ U 

°2  ~ + 

(K  ,**  + K+) n 


and  A = (K  **  + k.+) (K  **  + K *)/(K+K  **  + K *K+  + K**K  *) 


It  can  be  shown  that  Eq.  (30)  is  sufficient  to  ensure  the  validity 
of  Eq.  (29a-d) . In  such  cases  Eq.  (28)  would  predict  temporally 
decaying  solutions  indicating  stable  discharges.  The  parameter 
6^  and  62  are  the  ratios  of  characteristic  time  for  loss  of  elec- 
trons by  electron  impact  to  that  due  to  heavy  particle  collisions 
for  the  two  excited  states.  It  is  observed  that  when  §2  or  6^ 
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approaches  zero  (with  A6^  or  A6^  » 1),  two-step  ionization 
limit  similar^ 17 ^ to  that  derived  by  Daugherty,  Mangano  and 

(9) 

Jacob  is  recovered.  This  corresponds  to  a situation  when 
n^  and  n2  are  completely  mixed  by  electrons.  When  6^  and  <S2 
both  approach  zero,  the  usual  single-step  stability  criteria 
(6  > v.  , and  an  > v.  ) is  recovered.  For  this  case  the  ex- 

xo  e0  x0 

cited  state  population  will  come  into  thermal  equilibrium  with 

(18) 

the  secondary  electrons.  ' It  is  noticed  from  the  expression 
for  e that  when  electron  mixing  of  the  excited  states  is  unim- 
portant (i.e.,  6^,62  > 1)/  its  value  is  bounded  (note  that 
A > 1)  by  0 < e < 1.  The  stability  criteria  for  three-step 

dominant  ionization,  therefore,  is  given  by  an  (or  3)  > mv . 

e0  x0 

where  1 < m < 3/2  for  recombination  and  2 < m < 3 for  attach- 
ment dominated  discharge.  Our  numerical  calculations  for  laser 
discharges  typically  indicate  6^  ~ 2,  S2  1 . 0 and  A s 1.1  for 
CdHg*  and  6^  « 13,  S2^l  and  A - 2 for  HgCl*  discharge.  These 
correspond  to  values  of  m = 1.1  and  2.5,  respectively.  In 
general,  knowledge  of  the  parameters  6^,  6^,  A and  the  equilibrium 

values,  n , v.  , etc.  (obtained  by  solving  Eqs.  (18)  to  (20), 
e0  x0 

would  determine  the  value  of  m below  which  the  discharge  will  go 

unstable.  Note  also  that  when  the  electron  production  rate 

S ^ = 0,  then  from  Eq.  (18)  the  equilibrium  condition  for  an 

attachment  dominated  discharge  is  8 = v.  , which  is  in  violation 

x0 

of  Eq.  (30)  when  e > -1.  Hence,  an  external  source  for  electron 
production  is  essential  for  discharge  stability  of  an  attachment 
dominated  discharge  when  e > -1.  By  a similar  reasoning,  a 


] 


i 
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nonzero  S ^ is  essential  for  the  stability  of  a recombination 
dominated  discharge  when  e > 0. 

In  order  to  demonstrate  the  validity  of  the  linearized 
analysis  presented  earlier,  we  have  also  numerically  solved  the 
system  of  Eqs . (18)  to  (20).  Figure  5 shows  the  results  of  such 

an  analysis.  Discharge  conditions  were  selected  such  that 
m ^ 3 for  this  case.  On  the  left-hand  side  we  have  the  stable 
discharge  condition,  i.e.,  one-third  of  the  attachment  rate  is 
slightly  greater  than  the  equilibrium  ionization  rate.  Note 
that  the  discharge  current  reaches  a constant  value  asymptotically. 
If  we  keep  everything  constant  but  decrease  the  attachment  rate 
by  5%,  we  observe  that  the  ionization  rate  increases  and  after 
about  0.8  ysec  becomes  greater  than  one-third  of  the  attachment 
rate.  The  discharge  current  grows  and  shows  stronger  than  ex- 
ponential growth  after  about  0.8  ysec.  In  Figure  6 we  show  a 
similar  behavior  for  recombination  dominated  discharge.  In  this 
case,  discharge  conditions  were  selected  such  that  m 1.5.  On 
the  left-hand  side  we  have  again  the  stable  discharge  condition 
as  long  as  ang  >3/2  and  the  discharge  current  reaches  a con- 
stant value  asymptotically.  If  the  recombination  rate  is  reduced 
by  10%,  we  observe  that  the  ionization  rate  increases  and  becomes 
> 3/2  ang  at  about  0.7  ysec.  The  discharge  current  exhibits  a 
stronger  than  exponential  growth  after  about  0.7  ysec. 

The  volumetric  instability  described  in  this  letter  predicts 
a temporal  increase  of  the  secondary  electron  dens.’  Ly  that  is 
faster  than  exponential  in  a spatially  uniform  manner.  For  spatial 
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Figure  5 Results  of  a Numerical  Analysis  of  a System  of 
Eqs.  ( 1 ) — ( 3 ) for  Attachment  Dominated  Discharge 
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